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ABSTRACT 

We show that the circular character of continuum structures observed in the coma of comet 
Hale-Bopp around the perihelion passage is most likely due to a dust jet from a large extended 
active region on the surface. Coma morphology due to a wide jet is different from that due to a 
narrow jet. The latter shows foreshortening effects due to observing geometry, wider jet produces 
more circular features. This circularization effect provides a self-consistent explanation for the 
evolution of near-perihelion coma morphology. No changes in the direction of the rotational 
angular momentum vector are required during this period in contrast to the models of Schleicher 
et al. This circularization effect also enables us to produce near-circular coma features in the S-E 
quadrant during 1997 late February and therefore questions the basic premise on which Sekanina 
bases his morphological arguments for a gravitationally bound satellite nucleus. 

Subject headings: comets: general — comets: individual (Hale-Bopp C/1995 01) — methods: numerical 


1, Introduction 

A primary characteristic of the evolving coma 
morphology of comet Hale-Bopp (C/1995 01) was 
the presence of repetitive near-circular dust fea- 
tures around the time of the perihelion passage 
(e.g., Mueller et al. 1999, and references therein). 
These features which persisted for nearly a three 
month period centered around the perihelion pas- 
sage of 1997 April 1 showed a periodicity of « 11.3 
hrs. This periodicity can be identified with the 
rotation of the nucleus (Farnham et al. 1998; Li- 
candro et al. 1998; Jorda et al. 1999; Warell et 
al. 1999). 

During the near-perihelion months, images cor- 
responding to emission features CN, C 2 , C 3 , NH, 
and OH (Larson et al. 1997; Schleicher et al. 
1998; Lederer et al. 1999; Mueller et al. 1999) 
show complete spiral structures around the nu- 
cleus. Careful image enhancement by Schleicher 
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et al. (1998) has revealed that the morphology in 
the continuum images also shows spiral behavior. 
But the signature therein is only barely evident 
in the anti-solar direction. These spiral structures 
were nearly circular Archimedean spirals (i.e., dis- 
played little or no foreshortening effects). They 
maintained approximately the same shape during 
the three month period centered around the peri- 
helion when the Earth-comet-sun (E-C-S) geome- 
try changed significantly. 

Attempts have been made to determine whether 
the nucleus is in a principal axis state or not as 
well as to fix the direction of the rotational angu- 
lar momentum vector, M (Sekanina 1998a; Jorda 
et al. 1999; Licandro et al. 1999; Samarasinha 
et al. 1999; Vasundhara & Chakraborty 1999). 
Simulations aimed at replicating near-perihelion 
coma morphology have encountered many difficul- 
ties (Sekanina 1998b; Vasundhara & Chakraborty 
1999) and Sekanina has argued that two gravita- 
tionally bound active nuclei might explain the 
evolution of near-perihelion coma morphology. 

Schleicher et al. (1998, 1999), in their attempts 
to mimic the evolution of near-perihelion coma 
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morphology, deduced a greater than 90° preces- 
sion of the axis of the “cone” generated by the rel- 
evant active region on the rotating nucleus. This 
precession has occurred over a three month period. 
Furthermore, this “cone axis” , which can be iden- 
tified with the vector M in the case of a principal 
axis rotator, approximately followed the direction 
of Earth which also moved nearly 90° over the said 
period, an unlikely or extremely fortuitous situa- 
tion. 

In this Letter we present new coma simulations 
that include a hitherto unexplored projection ef- 
fect which gives a self-consistent explanation for 
the evolutionary behavior of near-perihelion coma 
morphology without involving exceptional physi- 
cal or dynamical situations. 

We first show (section 2) that the coma mor- 
phology due to observationally compatible com- 
plex spin states are qualitatively the same as that 
due to a principal axis rotator. In section 3 we de- 
termine the probable rate of forced precession of 
M due to outgassing from the nucleus and show 
that it is unlikely to have moved more than few 
degrees during the 1997 March-May period. In 
section 4 we calculate the total area of the ac- 
tive regions on the nucleus and demonstrate that 
the dominant active region must have considerable 
angular extension (« 40°). Numerical simulations 
that emulate the coma morphology due to such an 
active region are then presented that show that a 
precession of M is not required to explain what 
is observed. Section 5 discusses certain additional 
issues which are of relevance to coma morphology. 

2. Can a Non-Principal Axis Spin State 
Explain the Evolution of Coma Mor- 
phology? 

During the 1997 March-May period, the near- 
circular spiral pattern of the coma morphology 
was nearly unchanged whereas the Earth direction 
(and the solar direction) moved by nearly 90°. To 
explain this Schleicher et al. (1998, 1999) pro- 
posed a continual precession of the vector M over 
this three month period. 

An immediate question that arises is whether 
a non-principal axis rotation state for the nu- 
cleus of comet Hale-Bopp can provide a more sim- 
ple explanation for this behavior. Samarasinha 
et al. (1999) identified two categories of non- 


principal axis spin states which can reproduce reg- 
ular repetitive coma structures. They are (a) a 
spin state with a small precession angle of the or- 
der of few degrees, and (b) a spin state in which 
one component period is several orders of magni- 
tude larger than the other period. For this lat- 
ter case, the shorter period corresponds to the 
precession of the angular velocity vector around 
the angular momentum vector, M (Ames & Mur- 
naghan 1929; Samarasinha et al. 1999). For 
both of these complex spin scenarios, the angu- 
lar velocity vector, the angular momentum vec- 
tor, and the “cone axis” must be aligned within 
few degrees of each other. As in the case for the 
principal axis states, in all observationally com- 
patible non-principal axis states, the “cone axis” 
can for all practical purposes be approximately 
identified with the vector M. Therefore, if a pre- 
cession of the vector M is invoked for a princi- 
pal axis rotator to explain the evolutionary behav- 
ior of near-perihelion coma morphology , the same 
holds true for the observationally compatible com- 
plex spin states. 

This does not mean there are no morphologi- 
cal differences between coma structures generated 
by a principal axis rotation and an observation- 
ally compatible complex rotation. Only that the 
“cone axis” will in both cases be nearly aligned 
with M all the time. For example, the rotation as- 
sociated with the longer period of a complex spin 
state with a large ratio between component peri- 
ods (the second complex spin scenario described 
earlier), would result in variations in the opening 
angle of the “cone” with time. This may manifest 
itself as a time dependent (related to the longer 
component period) variations in the projected out- 
flow velocity of the coma structures. Thus, irre- 
spective of whether a simple or a complex rotator 
is involved in the type of modeling done by Schle- 
icher et al. (1998, 1999), the result will imply a 
large precession of M. 


3. Dynamical Constraints on the Evolu- 
tion of Angular Momentum 


The e- folding time scale, tm> for changes in M 
due to reaction torques caused by outgassing is 
given by (Samarasinha et al. 1986; Jewitt 1992) 


7m = 


dVm 


( 1 ) 
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where P(= 27r/cj) is the rotation period, p is the 
bulk density of the nucleus, r is the nuclear radius, 
d is the effective moment arm for the net torque, 
V is the gas outflow velocity at the nucleus, m is 
the gas production rate, and k is a constant with 
a value of about 2. 

If d is taken to be r, we obtain a strict lower 
limit to tm- By adopting P = 11.3 hr, p = 0.4 g 
cm -3 , r = 30 km (Weaver & Lamy 1999; Y. 
Fernandez et al., in preparation), V = 1 km 
s~\ and m — 3 x 10 8 g s“ x (corresponding to 
Qh 2 o = 10 31 molecules s -1 (Biver et al. 1999)), 
we arrive at tm > 10 yrs. A more realistic value 
for d may be smaller by about 1 order of magni- 
tude (Jewitt 1999) and therefore tm is more likely 
to be of the order of 100 yrs. 

In addition, prior studies on spin state evolution 
show that large changes in M, due to outgassing 
torques, are accompanied by correspondingly large 
changes in the rotation period P (Samarasinha & 
Belton 1995). However, in the case of comet Hale- 
Bopp, no such change in the 11.3 hr periodicity 
during the 1997 March-May period was detected 
(Schleicher et al. 1998). Therefore, it is most un- 
likely that the direction o/M would have changed 
by more than few degrees during the approximate 
three month period around the perihelion. 

4. Coma Morphology Due to an Extended 
Active Region 

4.1. Size of Active Regions on Hale-Bopp 

We now examine whether the existing coma 
simulations, all based on point source or narrow 
jet models, are reasonable. 

First, we estimate the size of a single active 
region which could be responsible for the domi- 
nant coma feature. We adopt a perihelion water 
production rate of 10 31 molecules s -1 (Biver et 
al. 1999) corresponding to a perihelion distance of 
0.91 AU. Based on simple thermal models (Cowan 
&: A’Hearn 1979; Julian et al. 1999), this corre- 
sponds to an area of approximately 2 x 10 3 km 2 for 
the active region, i.e., an active region (or group of 
regions) which occupies 10% or more of the total 
surface area if the upper limit to the nuclear ra- 
dius is taken to be 40 km (Weaver & Lamy 1999). 
If only a fraction, say 40%, of the dust production 
occurs in the single dominant active region, then 
it would have an angular extension of the order 


of 40°. The exact angular extensions would, of 
course, depend on the shape of the active region 
and the actual nucleus radius itself. 

This contrasts with the existing simulations of 
coma features in Hale-Bopp, including those by 
the author himself, which have usually been car- 
ried out assuming either a point source model for 
the active region or a narrow jet. The physical 
situation questions this approach and we reexam- 
ine the coma morphology using a large extended 
active region as the basis of the modeling. 

4.2. Simulations of Coma Morphology 

Figure 1 shows the dependence of coma mor- 
phology for point source dust jets with different 
gaussian angular cross sectional widths. For all 
these cases, the observer is inside the “cone” swept 
out by the jet axes due to the rotation of the nu- 
cleus, the scenario corresponding to the near per- 
ihelion time. Since we are simulating coma struc- 
tures with an age of the order of only several rota- 
tion periods, and we are interested in identifying 
the behavior of coma morphology as a function 
of the jet width, radiation pressure effects are not 
considered. It is found that the spiral structure be- 
comes more circularized as the jet width increases, 
i.e., as the jet width increases, the foreshortening 
effects are reduced and the perceived length of the 
“minor axis 77 * approaches that of the “ major axis 7 ’. 

The following qualitative explanation will fur- 
ther clarify this non-intuitive jet- width dependent 
projection effect. 

For a narrow jet, foreshortening occurs in the 
“short axis” direction defined by the intersection 
of the sky plane and the plane which contains the 
line-of-sight (LOS) and M. Consider dust grains 
ejected during a short interval when the jet is in 
the plane containing the LOS and M. The resul- 
tant image would show a small dust patch in the 
direction of the “short axis”. The apparent out- 
flow velocity for that dust patch would be V& sin#, 
where Vd is the radial dust outflow velocity and 
6 is the angle between the LOS and the direction 
at which dust was ejected. Now, consider a wide 
jet where all these ejected dust grains are spread 
over a shell of radius R defined by the dust outflow 
velocity Vd- If the angular distribution of dust is 
sufficiently wide, the peak brightness in an image 
in the direction of the “short axis” would occur 
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further out than in the narrow jet because of the 
extension of the shell along the LOS- This descrip- 
tion crudely illustrates how the foreshortening ef- 
fects are reduced when you have a wide jet. 

An alternative way of demonstrating this effect 
is shown in figure 2 where we show the coma mor- 
phology due to dust from an extended active re- 
gion. To simulate this morphology, we co-added 
five point source model jets which are placed equi 
distant within the extended active region. Images 
in figure 2 mimic the dust outflow from an ex- 
tended active region on Hale-Bopp at two widely 
separated times and show how a nearly circular 
coma morphology is maintained. Despite large 
changes in the E-C-S geometry, this morphology 
maintains its circular character over time scales of 
months around the perihelion time consistent with 
observations. 

This circularization of coma features by a wide 
jet provide a natural self-consistent explanation for 
the structures observed in Hale-Bopp over approxi- 
mately three months around the perihelion without 
requiring any precession of M. 

In addition, as we show in figure 3, outgassing 
from a wide extended area, unlike a narrow jet, 
generates the near-circular coma structures in the 
S-E quadrant during 1997 February 23 for the M 
direction assumed by Sekanina (1998b). This re- 
sult questions the basic premise upon which he has 
argued for two active nuclei for comet Hale-Bopp. 

5. Discussion 

5.1. Radial Outflow Velocity of Dust 

A comparison of coma features due to jets of 
different widths (figure 1) indicate that the circu- 
larization occurs as a reduction of foreshortening 
effects in the direction normal to the maximum 
width of the feature. For wide jets, the maxi- 
mum width of the features corresponds to the ra- 
dial outflow velocity of dust. Therefore, separa- 
tion between repetitive structures in the direction 
of maximum width is a direct measure of the radial 
dust outflow velocity. 

5.2. Other Factors Affecting Coma Mor- 
phology 

The simulations presented in this Letter as- 
sume that the dust jets maintain the same angular 


cross sectional widths irrespective of cometocen- 
tric distance. Hydrodynamical simulations, espe- 
cially that close to the nucleus (Kitamura 1990; 
Crifo & Rodionov 1999, and references therein) in- 
dicate the possibility of substructures due to shock 
formation if there are several spatially separated 
source regions. In addition, if the whole active re- 
gion is not outgassing uniformly but has localized 
high activity at certain parts, then that could also 
be responsible for some of the substructure. It is 
not clear how these effects will be propagated out 
to the observable coma and relevant calculations 
are yet to be carried out, a task beyond the scope 
of this Letter. Here, we assume that they would 
only cause secondary effects and indeed some of 
the fainter substructures seen in the coma mor- 
phology of comet Hale-Bopp may be related to 
such effects. A detailed examination of hydrody- 
namical effects close to the nucleus and their possi- 
ble influence on the observable coma morphology 
will be discussed elsewhere (M. Belton et al., in 
preparation). 

An increase in the angular cross sectional width 
of dust jets in the circumnuclear region as a 
function of increasing cometocentric distance is 
demonstrated by Combi et al. (1999) based on 
their hydrodynamical simulations. Therefore, the 
angular cross sectional width of a dust jet at ob- 
servable distances may in fact be larger than its 
original angular width at the nucleus, again justi- 
fying the choice of a wide jet. 

5.3. Coma Morphology Due to Multiple 
Active Regions 

In addition to the spiral coma structure due to 
a single dominant active region, there are at least 
one (Schleicher et al. 1998) or perhaps two (Li- 
candro et al. 1998; Farnham et al. 1999) jet 
features which intersect with the primary feature 
at shallow angles for each repetitive cycle. As nar- 
rower jets would show different curvatures than 
the curvature of the wider dominant jet, these 
secondary jets are likely to be originating from 
smaller active regions on the surface. These need 
to be taken into account during future detailed 
coma modeling. 
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6. Summary 

Coma morphology due to a wide jet is shown 
to be substantially different from that due to a 
narrow jet. The morphology due to a narrow jet 
shows strong foreshortening effects due to observ- 
ing geometry, while a wider jet produces more cir- 
cular features. 

This jet-width dependent projection effect has 
the following immediate implications. 

(1) The nearly unchanged repetitive circu- 
lar spiral structures observed over approximately 
three months in Hale-Bopp can be explained most 
simply by a wide jet originating from a large ex- 
tended active region on the nucleus. Little, or no, 
change in the direction of M during this time in- 
terval is required. This scenario is consistent with 
dynamical calculations. 

(2) Contrary to the claim of Sekanina (1998b), 
simulations indicate the existence of near-circular 
coma features in the S-E quadrant during 1997 
late February. This seriously questions the basic 
premise upon which he bases his morphological 
arguments for two active nuclei. 

(3) The separation between the repetitive cycles 
of the coma structure in the direction of maximum 
width is a direct measure of the radial dust outflow 
velocity. 

I thank Dr. M. Belton for many helpful sugges- 
tions, Dr. B. Mueller for going through an early 
draft, and Drs. D. Schleicher and Z. Sekanina for 
discussions on their work. I thank the referee, 
Dr. T. Farnham for his valuable comments. This 
work was supported by NASA Planetary Astron- 
omy Program. 
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Fig. 1. — Simulated coma morphology as seen 
from Earth for the Earth-comet-sun geometry of 
1997 March 13 for dust jets having standard de- 
viation, (j , of 1°, 5°, 15°, and 30° (clockwise from 
top left) for the gaussian cross sectional width. For 
all four cases, the other parameters were kept un- 
changed. Earth is inside the “cone” swept out by 
the axes of the jets. For each panel, the nucleus 
is at the center. North is up and east is to the 
left. A single dust outflow velocity is assumed. 
The dust production rate is assumed to be pro- 
portional to (0.1-fcosa) during day time and 0.1 
during the night, where a is the solar zenith angle. 
This night time emission, which helps the eye to 
trace the spiral structure, is primarily confined to 
an approximately 90° sector in the westerly direc- 
tion. 

Fig. 2. — Simulated coma morphology as seen 
from Earth for E-C-S geometries of 1997 March 13 
(left) and 1997 May 6 (right) for otherwise iden- 
tical parameters. Point source dust jets of a of 
15° originating from a five point equi-distant grid 
located within an extended active region are co- 
added in each image. The five points are located 
at the center of the active region (same physical 
location as that in figure 1), at longitude displace- 
ments of ±20°, and at latitude displacements of 
±20°. For both times, Earth is inside the “cone” 
swept out by the axis of the central jet. Same dust 
production rate as in figure 1 is used. Spatial cov- 
erage in May image is larger since the geocentric 
distance is larger. 

Fig. 3. — Simulated coma morphology as seen 
from Earth due to a jet of a of 1° (left) and that 
due to an extended active region (right) for 1997 
February 23 for the same M direction of Sekan- 
ina (1998b). We used a five point grid similar to 
that described in the caption for figure 2 to mimic 
the dust outflow from an extended active region. 
The origin of the point source jet on the left as 
well as the center point of the extended active re- 
gion are located at latitude + 10°. For each panel, 
the nucleus is at the center. North is up and east 
is to the left. Zero dust production is assumed 
during the night time. Notice that the image on 
the left corresponding to a narrow jet is much like 
what Sekanina (1998b, figure 3) obtained, while 
an extended active region produces near-circular 
features in the S-E quadrant. 
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